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Abstract 

We study the phenomenological consequences of amphtude-corrected post-Newtonian (PN) grav- 
itational waveforms, as opposed to the more commonly used restricted PN waveforms, for the 
quasi-circular, adiabatic inspiral of compact binary objects. In the case of initial detectors it has 
been shown that the use of amplitude-corrected waveforms for detection templates would lead 
\^ ', to significantly lower signal-to-noise ratios (SNRs) than those suggested by simulations based ex- 

I clusively on restricted waveforms. We further elucidate the origin of the effect by an in-depth 

' analytic treatment. The discussion is extended to advanced detectors, where new features emerge. 

^ ■ Non-restricted waveforms are linear combinations of harmonics in the orbital phase, and in the 

O . frequency domain the kth harmonic is cut off at kf^so, with fLSO the orbital frequency at the 

! last stable orbit. As a result, with non-restricted templates it is possible to achieve sizeable signal- 

lO I to-noise ratios in cases where the dominant harmonic (which is the one at twice the orbital phase) 

' does not enter the detector's bandwidth. This will have important repercussions on the detection 

^ ■ of binary inspirals involving intermediate-mass black holes. For sources at a distance of 100 Mpc, 

. taking into account the higher harmonics will double the mass reach of Advanced LIGO, and that 

, of EGO gets tripled. Conservative estimates indicate that the restricted waveforms underestimate 

I detection rates for intermediate mass binary inspirals by at least a factor of twenty. 
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I. INTRODUCTION 



A number of interferometric gravitational-wave detectors (LIGO, VIRGO, GEO600 and 
TAMA) are now running and the LIGO detector has reached design sensitivity. In 
the coming decade the Initial LIGO detector will be upgraded to Advanced LIGO. For all 
such detectors, the inspiral signals from neutron star and/or black hole binaries are among 
the strongest expected sources (see e.g. jl] for a review). In the late stages of inspiral, 
gravitational radiation backreaction will have circularized the components' orbits, and before 
the final plunge there will be an "adiabatic" regime where the period of a single orbit is much 
shorter than the inspiral timescale. This part of the inspiral is relatively "clean" and well- 
understood, so that the waveforms it produces can be modeled with great precision. This 
has been done in the so-called post-Newtonian (PN) approximation to general relativity, 
where waveforms are expressed as expansions in the orbital velocity v (see jsf for a review 
and extensive references). 

The best waveforms currently available are of order in amplitude 0] and v'^ in phase 
jsl, which in the usual notation corresponds to 2.5PN and 3.5PN orders, respectively. These 
take the form of linear combinations of harmonics in the orbital phase. The second harmonic 
is the one that will usually dominate; it is the only one with a zeroth-order PN contribution 
to its amplitude. In most of the literature on adiabatic inspiral, one uses the restricted 
post-Newtonian approximation j^, where all amplitude corrections are discarded and only 
PN contributions to the phase are taken into account. Thus, the restricted PN waveform 
consists of just the second harmonic with a prefactor at zeroth PN order. 

An efficient way of searching for inspiral signals in data is matched filtering 0], which 
involves a bank of templates. Almost all templates currently used are based on the restricted 
PN approximation. The same kinds of waveforms that go into template banks are also used 
as simulated signals injected into stretches of data to evaluate algorithms that search for 
real events and veto spurious ones. These include the standard restricted PN waveforms, 
the so-called P-approximants [sl, and the effective one-body (EOB) waveforms j^. The 
latter two result from resummation schemes designed to improve on the convergence of the 
phasing. The phenomenological templates for binary black hole inspiral that were recently 
proposed by Buonanno, Chen, Pan, and Vallisneri [lO|] do have some simple corrections to 
their amplitudes, but they are still quite limited compared to the full amplitude-corrected 
waveforms. Overall the emphasis has been on the evolution of the phase rather than the 
amplitude, due to the belief that it is more important to know the phasing and number 
of cycles (or more precisely the number of useful cycles [ll[) of the signal in the detector's 
bandwidth. 

To our knowledge it has never been investigated in detail to what extent the use of 
restricted waveforms as templates and simulated signals is justified beyond heuristic con- 
siderations. Sintes and Vecchio [l2| did study the effect of 0.5PN amplitude corrections on 



signal-to-noise ratios (SNRs). However, as shown in [13|, this happens to be a very special 
case; in going from OPN to 0.5PN in amplitude, SNRs will increase, but from IPN onwards 
there will be a pronounced drop, at least for stellar mass inspirals as seen in initial detectors. 
The effect of 2PN amplitude corrections on parameter estimation was studied by Hellings 
and Moore [13], but specifically in the context of LISA and without discussing the effect 
of higher harmonics on detection rates for systems with various masses. Here we inves- 
tigate the use of non-restricted versus restricted waveforms for the purposes of detection, 
focusing on three ground-based detectors: Initial LIGO, Advanced LIGO, and a possible 
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third-generation detector called the European Gravitational- Wave Observatory (EGO) |l5| . 

Given a waveform h, the signal-to-noise ratio (SNR) for the "detection" of h using its 
normalized counterpart as a template is given by p[h] = (/i|/i)^/^, where ( . | . ) is the usual 
inner product in terms of the noise power spectral density Sh{f) of the detector. With the 
convention x(/) = J^^x(t) exp{—2'JTift)dt for Fourier transforms, 

Let ho be a restricted waveform and h a waveform that is of 2.5PN order in amplitude 
and 3.5PN in phase, with both having the same parameters. With the above notation, 
P[ho] is the kind of SNR one encounters in current simulated searches. In the future one 
may want to consider amplitude-corrected waveforms h both for templates and simulated 



signals, in which case one would arrive at SNRs p[h]. As indicated in an earlier paper [13 
in the context of initial detectors and stellar mass binaries, one has 



p[h]<p[ho]. (1.2) 

Hence, surprisingly, if one were to use the best available amplitude-corrected waveforms 
for detection templates, one should expect SNRs in actual searches to be lower than those 
suggested by simulations based purely on restricted PN waveforms. In Initial LIGO this 
overestimation can be as large as 25% depending on what the parameter values are. Because 
SNR is inversely proportional to distance, this corresponds to an overestimation of the 
accessible spatial volume, and hence the detection rate, by up to a factor of two. Moreover, 
SNRs exhibit a downward trend as the post-Newtonian order of the amplitudes is increased 
from IPN to 2.5PN in steps of 0.5PN. Should this trend continue then the overestimation 
would be worse still when taking into account amplitude corrections beyond the highest PN 



order currently available. A rough indication of the origin of (11.21) was already given in 13 
The effect can be traced to the PN amplitude corrections of the dominant harmonic. Here 
we will provide a much more detailed analysis. 

The situation is different for second and third generation detectors. In the frequency 
domain, the third and higher harmonics in amplitude- corrected waveforms are cut off at 
higher frequencies than the dominant one. Because the cut-off frequencies are inversely 
proportional to total mass, this means that higher harmonics may still enter the detector's 
bandwidth even if the dominant harmonic does not. In advanced detectors these harmonics 
can lead to a sizeable SNR. If the second harmonic does not enter the bandwidth, or if it 
enters only in a small frequency interval, the inequality (11.21) is reversed. This opens up the 
possibility of seeing intermediate-mass inspiral events with higher total mass than one would 
expect on the basis of the restricted waveforms. Indeed, the inclusion of amplitude-corrected 
higher harmonics in detection templates would effectively double the accessible mass range 
of Advanced LIGO, and it would triple that of EGO. For intermediate-mass inspirals the 
detection rates from non-restricted waveforms are several orders of magnitude larger than 
the ones computed from the restricted waveforms. Stellar mass inspirals would be detectable 
with EGO throughout a significant part of the visible Universe. 

This paper is structured as follows. In section [TTl we discuss the amplitude-corrected wave- 
forms, as seen in a detector, up to 2.5PN in amplitude and 3.5PN in phase in the stationary 
phase approximation. In section IIIII we compare signal-to-noise ratios using restricted and 
non- restricted waveforms as templates and simulated signals in Initial LIGO. The numerical 
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results are explained analytically up to IPN in amplitude. We then focus on the advantages 
of non- restricted templates in advanced detectors (section HVll . including a third generation 
detector called EGO. Finally we summarize our results and conclusions. Expressions for the 
harmonics in the amplitude-corrected waveform can be found in Appendix A. We will go 
into some detail, outlining the general form of the harmonics up to 2.5PN in amplitude and 
givin g ex plicit expressions up to IPN.^ These expressions are also to be the basis for future 
work [16|, where we will use the amplitude-corrected waveforms for parameter estimation. 



Sensitivity curves for detectors will be taken as in [17|, except for the EGO noise curve 
which can be found in Appendix C. Unless stated otherwise we use units such that G = c = 1. 
We use the notation (p, g)PN for waveforms that have PN corrections up to order p in 
amplitude and q in phase. Integrals of the form (11.11) can not be performed analytically and 
were obtained numerically using the software package Mathematica. Finally, when plotting 
quantities against mass in a cosmological context, we will give precedence to physical mass 
over (redshifted) observed mass. This is because we will encounter relatively large redshifts, 
in which case the use of observed mass would obscure the astrophysical consequences of the 
results. 



II. AMPLITUDE-CORRECTED POST-NEWTONIAN WAVEFORMS 

The waveforms in the two polarizations take the general form 

/^+,x = {<i + xV2</f ^ + xH^^l + ^ + x'H^Z + ^'^'^x ^ } (2-1) 

where r is the distance to the binary, M its total mass, and t] the ratio of reduced mass to 
total mass. The post-Newtonian expansion parameter is defined as x = = w^, 

with F{t) the instantaneous orbital frequency. The coefficients H^^^\ n = 0,...,5, are 
linear combinations of various harmonics with prefactors that depend on the inclination 
angle l of the angular momentum of the binary with respect to the line of sight as well as on 
rj] their explicit expressions can be found in jj]. The measured signal also depends on the 
polarization angle and the position in the sky through the detector's beam pattern functions 

h{t) = F^h+{t) + F^h^it). (2.2) 

Note that for ground-based detectors, which are the ones we will be concerned with, it is 
reasonable to approximate -F+,x as being constant over the duration of the signal. They 
depend on angles {6,(l),ip), where (6', 0) determine sky position while ip is the polarization 
angle. The signal (12.21) is a linear combination of harmonics of the orbital phase \l/(t) with 

offsets (f{k,m/2y- 

Np 2p 

^(^) = Ak,m/2){t) COs{k^{t) + iPik,m/2)), (2.3) 

k=l m=0 

where the coefficients A(^k,m/2) are functions of (r, M, r], 6, 0, tp, l) multiplied by x^™"*"^)/^. The 
orbital phase \l/(t) is a series in x, which in the case of non-spinning binaries is known to 



The general expression for the complete waveform in the stationary phase approximation up to 2.5PN in 
amplitude is available upon request as a Mathematica notebook. 
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3.5PN order. The number of harmonics Np depends on the PN order in amphtude, p; at 
2.5PN one has Np = 7. 

More exphcitly, the waveform h{t) is found as follows: 



Substitute the expressions (12. ip into (12.21) . taking the prefactors H^I^\ri,L) to be as 
in and beam pattern functions 

F+{e, 0, ^) = ^ (1 + cos^(^)) cos(20) cos(2^) - cos(^) sin(20) sin(2?/^), 

Fx {e, 0, ^) = i (1 + cos2(^)) cos(20) sin(2^) + cos(^) sin(20) cos(2V^). (2.4) 

In the resulting expression for h{t), collect cosines and sines of multiples of the orbital 
phase to arrive at 

Np 2p 

^^^) = £ £ [«(fc,m/2)(^) cos(A;^(t)) + Pik,m/2){t) sin(A;^(t))] , (2.5) 

fc=l m=0 

where the (y(k,s), P{k,s) depend on (r, M, r], l, 9, 0, ip), as well as on time through x. 
Combine corresponding sines and cosines into simple cosines to arrive at (12.31) . with 



A(^k,s) = sign(a(fc,,)) ^ a^^^^^ + /?(\^) (2.6) 

and 



a 



During the inspiral phase one has \dlnA(^k,s)/dt\ <^ 1 and \kd'^'^ /dt'^\ ^ {kd^ IdtY, in 
which case one can use the well-known stationary phase approximation (SPA) [18|] for the 
Fourier transform of (12. Sp : 



5 . ....... , ^ , / /I 



^ Em=0^(fcW2) (^(j/)) 



I ( l-Kjt [-j\-k^\t\-j\]- <^(fc,^/2) - 7r/4 



exp 



z ( 27r/t, - 7r/4 + fc^ ( ^/ 



(2.8) 



A dot denotes derivation w.r.t. time and is the coalescence time. The function Fii) is 
the instantaneous frequency associated with the orbital phase, lixFii) = \&(t), and t(/) is 
defined implicitly by F{t{f)) = f. The expression for the "frequency sweep" F in terms of 
F is given in Appendix A. To 3.5PN order, the function ip{f) takes the form 
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where tpc is the orbital phase at coalescence, and for the coefficients ^pi we again refer to 
Appendix A.^ 

The SPA for the {p, 3.5)PN waveform is then 



hsPAif) = 



E 



exp 



27r/te - vr/4 + k^j ( -f 



(2.10) 

where [ . ]p denotes consistent truncation to pth post-Newtonian order (i.e., the "Newtonian" 
prefactor /"^^^ is taken outside and the remaining expression is expanded in (27rM/)^/^ up 
to (27rM/)2p/3). 

The SPA is just one method of approximating the Fourier transform; the main alternative 
is the fast Fourier transform (FFT). However, as shown in 11, 2^], the difference between the 
two in terms of overlap and SNRs tends to be only a few percent for restricted waveforms, at 
least in initial detectors, and the situation improves dramatically as the detector's sensitivity 
at low frequencies is increased. We also note that the conditions for the applicability of the 
stationary phase approximation on the harmonics become more favorable with increasing k. 
The reason why we prefer to work with the SPA is that it will give us direct analytic insight 
into the effect on detection rates of the presence or absence of amplitude corrections. 

Finally, the expressions for the respective harmonics should not be used up to arbitrarily 
large frequencies; some maximum orbital frequency must be chosen at which waveforms are 
truncated. In the limit where one of the component masses is sent to zero while keeping 
total mass fixed, there is a well-defined last stable orbit with orbital frequency 



LSO 



63/227rM' 



(2.11) 



where M is total mass. For binary systems with comparable component masses the situation 



is more involved [21[ , but for simplicity we adopt (12. lip as a formal expression also in that 
case. In the time domain we will consider the waveform to be valid only up to a time 
determined by F(t) = fiso- In the frequency domain this roughly corresponds to cutting 
off the kth harmonic at a frequency / = kf^so- To reflect this restriction on the validity 
of the harmonics, in practice we multiply the kth harmonic by 9{kfLso ~ f)i where 9{x) 
is the usual Heaviside function {9{x) = for x < and 1 otherwise). In the definition of 
the inner product ( . | . ), Eq. (11.11) . we then have fend = Tfiso, the frequency reach of our 
highest harmonic. The fact that higher harmonics have a higher frequency reach will have 
important implications, especially for advanced detectors. The lower cut-off frequency fs 
depends on the detector. For Initial LIGO, Advanced LIGO, and EGO, we will take it to 
be 40 Hz, 20 Hz, and 10 Hz, respectively (see Appendix C). 



^ Note that we are using slightly different conventions from e.g. [l^, where the focus was on the 
dominant harmonic and F{t) referred to twice the orbital frequency. Here primacy is given to the orbital 
quantities themselves, so that the dependence of F on and the frequency dependence of ip look a little 
different from what one often finds in the literature. 
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III. SIGNAL-TO-NOISE RATIOS: INITIAL DETECTORS 



We now turn to the issue of signal-to-noise ratio in Initial LIGO. Apart from binary 
neutron star inspirals, the inspiral events of main interest in the case of initial detectors 
are those involving stellar mass black holes. These result from massive main-sequence bi- 
nary stellar systems (field binaries) or from capture processes in galactic centers or globular 
clusters (capture binaries). Field binaries are expected to have total masses in the range 



10 Mq ^ M < 50 Mq while capture binaries may be somewhat heavier 22|]. Coalescence 



rate estimates for field binaries range from 
-3 in-i,.^-i 



10 * yr ^ to ~ 10 ^ yr ^ in our galaxy, corre- 



sponding to ~ 10 "^-10 ^yr ^ within a distance of 100 Mpc 0, 23, 24 1. For capture binaries 
lobular clusters, coalescence rates have been estimated at ~ 3 yr~^ within 600 Mpc 



m 

if 



26 



Below the implications of using non-restricted versus restricted waveforms for templates 
and simulated signals are first investigated numerically and then explained analytically. 



A. Numerical observations 



First let us look at the SNRs for inspirals at 20 Mpc in Initial LIGO. We consider three 
kinds of systems: a binary neutron star (NS-NS), a neutron star and a black hole (NS- 
BH), and a binary black hole (BH-BH). For concreteness we take the neutron star mass to 
be Mjvs = 1.4 Mq while for the black hole Mbh = 10 Mq. Table [T] displays the behavior 
of p[ho] and p[h] for {p, 3.5)PN waveforms h with increasing p, with Hq the corresponding 
restricted waveform. In going from p = to p = 0.5, for the asymmetric system p[h] 
increases slightly (as has been noted by Sintes and Vecchio [lH) while for the NS-NS and 
BH-BH there is no appreciable difference. However, at p = 1, p[h] drops below p[ho], and as 
p is further increased a clear downward trend is seen. At p = 2.5 one has p[h] < p[ho], with 
the two differing by 3.6%, 13.1%, and 21.0% for the NS-NS, NS-BH, and BH-BH inspirals, 
respectively. 





NS-NS 


NS-BH 


BH BH 


p 




p[h] 


p[ho] 


p[h] 


p[ho] 


p[h] 





6.465 


6.465 


13.492 


13.492 


30.928 


30.928 


0.5 




6.465 




13.936 




30.928 


1 




6.286 




12.563 




28.140 


1.5 




6.286 




12.417 




28.140 


2 




6.249 




12.091 




26.377 


2.5 




6.238 




11.933 




25.562 



Table I: Change in signal-to-noise ratios with increasing p in (p, 3.5)PN waveforms, for three dif- 
ferent systems at 20 Mpc. (Angles were chosen arbitrarily as 9 = (j) = it/6, ip = vr/4, l = '/r/3.) 



Thus, modeling templates and signals as restricted waveforms overestimates the SNRs. 
The fractional overestimation is given by 



p[h] 



-1, 



(3.1) 



where from now on we let h denote a (2.5, 3.5)PN waveform. 
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e 


(0,0) 


(0,0) 


0.129 


(vr/6,7r/3) 


(vr/6,0) 


0.127 


)5 


(vr/3,vr/8) 


0.132 




(vr/2,vr/4) 


0.141 


(vr/3,7r/3) 


(vr/6,0) 


0.127 


77 


(7r/3,7r/8) 


0.132 




(7r/2,7r/4) 


0.141 


(vr/2,0) 


(vr/6,0) 


0.127 


77 


(vr/3,vr/8) 


0.128 



Table II: Angular dependence of the ratio e for a NS-BH system as in Table HI 




Figure 1: A plot of e, the fractional overestimation of SNR due to the use of restricted PN wave- 
forms, as a function of rj and chirp mass. (The sudden drop in the bottom left corner is due to the 
fact that we imposed M > IMq.) In an important part of parameter space the overestimation is 
more than 20%. However, also note the negative values for very massive, asymmetric systems. 

In Table |T] specific choices were made for the angles [6, (f), l, ip) that encode the position 
and orientation of the source. Although p[ho] and p[h] will vary significantly as functions of 
these angles, this is not the case for e. Indeed, as can be seen in Table HTl for a NS-BH system 
the relative difference between p[ho] and p[h] has only a very weak angular dependence. 

In Fig. [1] we have plotted e as a function of i] and chirp mass M = Mri^/^. In a significant 
part of parameter space the overestimation is above 20%. Indeed, on the line corresponding 
to e = 0.2, lighter and heavier component masses mi, m2 are constrained as 1 Mq < mi < 
12 Mq and 12 < ma < 20 Mq. One has e < 0.25 for M < 30 Mq. The implied range 
of binary systems is of clear astrophysical interest. The overestimation shows a strong 
dependence on both the mass and the asymmetry of systems. 

We note that for systems that are both very massive and very asymmetric, e can become 
negative, and p[h] > p[hQ]. This occurs when the dominant harmonic enters the detector's 
bandwidth in only a small frequency interval. In that case both p[h] and p[ho] will be small, 
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at least for initial detectors; as we will show later on, this is not necessarily true for advanced 
detectors. 



B. Analytic treatment 



We have seen that in initial detectors, modeling signals as restricted waveforms generally 
leads to an overestimation of signal-to-noise ratio: 



p[h] < p[/io]. 



(3.2) 



We now analytically explain (13. 2p at IPN order in amplitude and within its domain of 
validity. 

The inequality (13. 2p is equivalent to {h\h) < (/io|/io)- Write h = Ylk=i ^'■'^^ "where the h'^^^ 
are amplitude-corrected harmonics. Leaving out the step functions in frequency, the latter 
are of the form 



J2Cik,m/2){2nMfy^/' 



m=0 



exp 



t [ 271 ft, - 7r/4 + kij[-f 



1 



(3.3) 



where C is a real function of chirp mass and distance while the coefficients C(fc,m/2) are 
complex functions of {9, L,r]). (Setting C(fc,m/2) = for m > would lead to the 
restricted waveform ho.) Now consider waveforms h for which /l5o ^ Is- In that case one 
can write 



{h\h) ^ ^(/i^'^^l/i^^)). 



(3.4) 



k=l 



because the various harmonics tend to interfere destructively with each other. (It can be 
checked numerically that the approximation is valid to within a few percent for stellar mass 
systems as seen in Initial LIGO.) The contributions (/I'^'^'^l/i'^'^^) as well as {ho\ho) take the 
form of integrals, the integrands of which all have a common factor 1/Sh{f)- For brevity, 
let h be of the (1,3.5)PN type, so that only the first four harmonics are present. Then up 
to order (27rM/)^/^, the remaining functions in the integrand can be written as 

Re[(^«(/))*/i«(/)] + Re[Ch^'\f)yh^'\f)] + Re[Ch^'\f)rh^'\f)] + Re[Ch^'\f)yU'\f)] 

= C^f -7/3 
X [K(2,0)|' 

+ (27rM/)2/3 (|C(i,i/2)r + IC(3,i/2)P + 2Re[C(2,o)C(2,i)]) 

+ (27rM/)^/3(|C(2,i)r + |C(4,i)r])]. (3.5) 

Now, 

|C(i,i/2)r + lC(3,i/2)P + 2Re[C(\o)C(2,i)] < (3.6) 

for any parameter values, and so one should expect {h\h) < {ho\ho), whence p[h] < p[ho]. 
For an explicit proof we refer to Appendix B. 

Obviously the inequality (13. 2p will not hold for systems where 2fLso < fs < Tfiso, 
because then (/lol^o) = while {h\h) > because of the contribution of the higher harmonics; 
recall that we have disregarded the step functions 0{kfLso~ f)- Also, if the second harmonic 
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enters the bandwidth in too small a frequency interval then the approximation (13.41) breaks 
down and the inequality f l3.2p may again be reversed; hence the regions of negative e in Fig. [H 
Both of these points will be important in the discussion of advanced detectors (section IIVI 
below). Finally, if h had been a (0.5, 3.5)PN waveform then in the inequality (13. 6p only the 
strictly positive contributions would have survived, so that once more p[h] > p[/io]- This is 
as observed by Sintes and Vecchio in [12] , and it is also evident in Table [H 



IV. ADVANCED DETECTORS 



We now turn to Advanced LIGO and EGO; for the former we use the estimated power 
spectral density from l3] while an analytic fit of the EGO PSD can be found in Appendix C. 
In addition to stellar mass binary inspirals, these will be able to see inspirals of intermediate- 
mass black hole binaries with total masses between 50 and more than a thousand solar 
masses. The latter may be formed in galactic nuclei in the process that leads to the formation 



of a supermassive black hole or in globular clusters (22|]. Coalescence rates are uncertain; 
assuming they are only 10~^ that of binary neutron stars 2^, there will be a few events per 
year within a distance of 2 Gpc. 

Below we first briefly discuss the effect of higher harmonics on the mass reach. Next we 
look at how they affect the redshift to which inspirals can be detected, after which we com- 
pute some estimates of detection rates with restricted and amplitude-corrected waveforms. 



A. Mass reach 




Figure 2: Plots of p\h] and /5[/io] as functions of total mass for Advanced LIGO (left panel) and 
EGO (right). Distance is fixed at 100 Mpc, and we assume mi/m2 = 0.1. Angles are as in 
Table HI At low masses one has 2fiso ^ fs and /o[/io] dominates. For sufficiently high masses, 
"^Ilso ^ fs, so that the dominant harmonic no longer enters the detector's bandwidth and the 
SNR for the restricted waveform vanishes. For such masses, higher harmonics in the amplitude- 
corrected waveform will continue to enter the bandwidth and can lead to significant SNRs. As a 
result, at the given distance the use of amplitude-corrected waveforms approximately doubles the 
mass reach of Advanced LIGO and triples that of EGO. 
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Here the emphasis will be on inspiral of intermediate mass binaries. It is generally 
thought that in a given cluster, usually only one intermediate mass black hole (IMBH) will 
form, which is left to interact with stellar mass compact objects j22l . 27| . Inspirals involving 
an IMBH will then tend to be asymmetric, with mi/m2 roughly in the range 10~^-10~^. 
In Fig. [5] we plot the SNRs p[h] and p[/io] against total mass M, keeping mi/m2 fixed at 
0.1 and taking the distance to be 100 Mpc. Here too, for low masses, modeling signals 
and templates as restricted waveforms would lead to a consistent overestimation of SNR: 
p[ho] > p[h]. However, as mass is increased, p[h] gradually comes to dominate. Eventually 
the k = 2 harmonic no longer reaches the detector's bandwidth, at which point p[ho] drops 
to zero. This occurs when 2fLso = (6^/^7rM)~^ = fs, with fs the detector-dependent lower 
cut-off frequency. However, for higher masses, the higher harmonics will still be in the 
bandwidth, and as can be seen from the plots of p[h], they can lead to a sizeable SNR. As 
a result, at a distance of 100 Mpc the use of amplitude-corrected templates approximately 
doubles the mass range accessible to Advanced LIGO, and it triples that of EGO! 

Let us assume event rates for intermediate mass inspiral to be in the order of a few per 
year within a distance of 2 Gpc. As can be inferred from Fig. [2] (by rescaling SNRs by 
a factor of 20), at that distance the SNR p[h] in Advanced LIGO will stay below 4 for 
mi/m2 = 0.1. By contrast, in EGO one has p[h] > 5 up to M ~ 600 Mq with a maximum 
of p[h] ~ 20 for M ~ 100 Mq, the approximate lower limit for intermediate-mass binaries. 
Although such inspirals will be accessible to Advanced LIGO, they would have to happen at 
distances of at most a few hundred Mpc to be seen with appreciable SNR, and coalescence 
rates are likely to be low within the corresponding spatial volume. This is much less of a 
problem for EGO, which has a significantly larger spatial coverage. 



B. Redshift reach 

In Fig. [21 distance was fixed at 100 Mpc, leading to large SNRs even when the dominant 
harmonic does not enter the bandwidth. We may now ask to what distances sources can 
be seen depending on how many harmonics reach the bandwidth. Since large distances will 
be involved, a proper treatment will have to take into account the effects of cosmological 
redshift z. This is achieved by making the following replacements in waveforms: 



r 



D{z), M — ^M' = {l + z)M, 



t,-^t',= il + z%, f^f = -J—^ (4.1) 
where in a flat Universe the luminosity distance D depends on redshift as 

with Hq the Hubble parameter at the current epoch while ^matter- and are dimensionless 
quantities defined on the basis of the total matter density and the cosmological constant, 
respectively. Wherever specific values are needed we will set Hq = 71 km/(Mpc s) and 
^matter = 0.27, Q/y = 0.73 Notc that the redshifting of frequencies will also affect 

the observed frequency of last stable orbit: /lso — fLso/{^ + z)- In general, the observed 
quantities are the redshifted ones. Nevertheless, when specifying systems we will continue to 
do so in terms of physical mass. In older work, precedence has usually been given to observed 
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mass. In the present study we will be dealing with relatively large redshifts, and in order 
to gain a clear understanding of the astrophysical implications of the results (in particular, 
we will want to distinguish between stellar mass and intermediate mass systems) it is then 
more convenient to work with physical mass. As an aside we note that assuming knowledge 
of Hq, ^matter ^ud VLk, and given a network of (at least three) detectors, luminosity distance 
can be measured @,Hl3l| and z can be solved for by inverting the relation (14. 2 p . Together 
with the observed mass, this leads to a value for the physical mass. 




Figure 3: The redshift reach of Advanced LIGO (left) and EGO (right) as functions of (physical) 
total mass for fixed SNRs of 10 with amplitude-corrected and restricted waveforms. We have fixed 
nil/ 1712 = 0.1, and angles are as in Table HI 



(OA) 




z (AdvLIGO, ho) 


z (AdvLIGO, h) 


z (EGO, ho) 


z (EGO, h) 


(0,0) 


(0,0) 


0.248 


0.198 


2.306 


2.157 


(vr/6,vr/3) 


(vr/6,0) 


0.191 


0.160 


1.865 


1.752 




(vr/3,7r/8) 


0.139 


0.130 


1.360 


1.320 




(vr/2,7r/4) 


0.098 


0.097 


0.942 


0.898 


(vrAvr/3) 


(vr/6,0) 


0.121 


0.100 


1.176 


1.057 




(7r/3,7r/8) 


0.087 


0.083 


0.829 


0.766 


» 


(V2,vr/4) 


0.058 


0.057 


0.524 


0.471 


(7r/2,0) 


(vr/6,0) 


0.114 


0.095 


1.107 


0.992 




(7r/3,7r/8) 


0.075 


0.071 


0.699 


0.641 



Table III: Redshift reach given a (9, 90) Mq system in Advanced LIGO and EGO with restricted 
waveforms ho and amplitude-corrected waveforms h, for p[ho] = 10 and p[h] = 10. We have again 
fixed mi/m2 = 0.1. 

In Fig. [3] we display the redshift at which inspirals can be seen as a function of total 
mass, still keeping mi/m2 = 0.1, with a fixed SNR of 10 for the restricted and amplitude- 
corrected waveforms, in Advanced LIGO and EGO. In the absence of an analytic expression 
for redshift reach as a function of SNR, redshifts were obtained using an elementary bisection 
method, producing successively better guesses for z until the value of the SNR differed from 
the desired one by less than 0.01%. In the Figure the angles are fixed; angular dependence 
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mi/ 1712 


z (AdvLIGO, ho) 


z (AdvLIGO, h) 


z (EGO, ho) 


z (EGO, h) 


0.2 


0.170 


0.154 


1.668 


1.608 


0.4 


0.204 


0.169 


1.965 


1.857 


0.6 


0.217 


0.175 


2.073 


1.942 


0.8 


0.222 


0.178 


2.113 


1.967 


1 


0.224 


0.179 


2.123 


1.975 



Table IV: Redshift reach for different mass ratios given a total mass of 100 Mq, in Advanced LIGO 
and EGO with restricted waveforms /iq and amplitude-corrected waveforms /i, for p[/io] = 10 and 
p[h\ = 10. Angles are again as in Fig. [3j 

is explored in Table UTTl Finally, Table HVl explores the dependence of the redshift reach on 
the ratio of component masses mi/m2. 

First consider Advanced LIGO. Because p[ho] tends to dominate over p[h] at low total 
mass (see Fig. [2]), fixing both to be equal to 10 leads to a larger redshift reach for the 
restricted waveforms at the low mass end. However, the redshift reach for the amplitude- 
corrected waveforms quickly catches up; after having peaked around M ~ 50 Mq it decreases 
relatively slowly. A dip is seen just before M ~ 220 Mq, where the dominant harmonic leaves 
the bandwidth; for the restricted waveform z drops to zero near that mass. From Fig. [3l 
Table llllt and Table IIVI we may conclude that with restricted and non-restricted waveforms 
alike, a large fraction of inspirals with M < 100 Mq will be detectable at redshifts of z ~ 0.05, 
and several times farther away depending on mass ratio, sky position, and orientation. 

With EGO the situation is qualitatively similar, but there the redshift reach is typically 
an order of magnitude higher than for Advanced LIGO. From Fig. [3], Table UTTl and Table ITVl 
we infer that a large fraction of inspirals with M < 100 Mq are detectable with p[/io] = 10 at 
redshifts z > 0.5, although this number can again be several times larger with better mass 
ratio, sky position, and/or orientation. Well positioned and oriented stellar mass inspirals 
are visible throughout a significant part of the visible Universe. 



C. Detection rates 

Phinney j32| estimated the number of stellar mass inspiral events out to a given distance 

-Kmsp — — -KmVK, l^.dj 

LmW 

where L{dmax) is the accumulated blue band luminosity for distances d < dmax, Lmw the 
blue luminosity of the Milky Way, and Rmw the coalescence rate in a Milky Way equivalent 
galaxy. Tutukov and Yungelson [1^ arrived at the estimate Rmw — 1-4 x 10~^yr~^ for 
binary black hole inspirals. From the L yon -Meudon extragalactic database j33| or the TuUy 



Nearby Galaxy Catalog [3J] one infers [35 



Ltoti^dfnax^ 1 n— 2 / dmax \ 



J^MW 

out to distances of a few hundred Mpc. 



10-^ j-™^ ^4_4N 

ViMpcy ^ ' 
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We are interested in obtaining detection rates over cosmological distances. For simplicity 
we will assume that the event rate per comoving volume is the same everywhere in the 
Universe. 

Consider black hole binaries at redshift z, inspiraling at a rate N = 3 /{An) x 
10~^ Mpc'^ Rmw per comoving volume per unit of cosmic time local to the event. As seen 



from Earth, the event rate between redshifts z and z + dz is 36 



where the comoving volume element per unit redshift is 

^^44 , («) 

with rc{z) the comoving distance; in terms of luminosity distance one has r^iz) = D{z)/{1 + 
z). The total event rate, as seen from Earth, up to redshift Zmax is then 







47riV r™"" D^{z)dz 



Ho Jo {1 + z)3^Q^atteril + + Qa'^ 



(4.7) 



where D{z) is given by (14. 2p . 

In the previous subsection we saw that stellar mass binaries will be detectable in Advanced 
LIGO at redshifts of z ~ 0.05 and much beyond. Setting Zmax = 0.05, we arrive at a 
(conservative) detection rate of 

Ri1:tr^0.ljT~\ (4.8) 

We have seen that EGO would be able to detect most stellar mass binaries with redshifts 
z > 0.5; setting Zmax equal to this value we arrive at a detection rate 

Rftei?ar^70yT-\ (4.9) 

corresponding to about one detection every few days. 

Following Fig. [3], for intermediate mass binaries with mass 50 Mq ^ M < 400 M©, the 
amplitude-corrected waveform leads to a redshift reach z > 0.3, while with the restricted 
waveform one has z > 0.1. Let RfMjuU ^e the detection rate from the amplitude-corrected 
waveform and RfM%str the one from the restricted waveform. If in each case we set Zmax 
equal to above values, 

- 20. (4.10) 

^IM,restr 

Needless to say, the above estimates are only indicative (and most likely overly conserva- 
tive); one could perform a more in-depth analysis along the lines of j35l . 33], but that would 



be outside the scope of this paper. Nevertheless, they already provide dramatic demonstra- 
tions of how important it will be to use the amplitude-corrected rather than the restricted 
waveforms for advanced detectors. 
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V. CONCLUSIONS AND FUTURE DIRECTIONS 



We have investigated the effect on signal-to-noise ratio and distance reach of the presence 
or absence of amphtude corrections in post-Newtonian gravitational waveforms for the quasi- 
circular, adiabatic inspiral of binary compact objects. This was studied using the stationary 
phase approximation, which allowed us to gain direct analytic insight. Our conclusions may 
be summarized as follows. 



For low-mass systems (M < 30 M©), which are the ones of main interest for initial 
interferometric detectors, the use of restricted waveforms as templates and simulated 
signals can lead to significantly larger signal-to-noise ratios and detection rates com- 
pared to what is attained by the best available, (2.5, 3.5)PN amplitude-corrected wave- 



forms. This runs counter to what one might naively conclude from earlier work [12 
where only 0.5PN amplitude corrections were taken into account. Indeed, in going 
from OPN to 0.5PN in amplitude, SNRs will increase; however, from IPN onwards we 
observed a steady decrease in the SNRs p[h] computed from the amplitude-corrected 
waveforms compared to those of the restricted ones, p[h^. The effect is the most pro- 
nounced for extreme to moderately asymmetric systems with large masses. For total 
mass up to 30 Mq the difference between p[h] and p[/io] can be as large as 25%, which, 
in initial detectors, corresponds to a difference of almost a factor of two in detection 
rates. We studied this effect analytically at the IPN level and found that, for the 
specified systems, it results from the amplitude corrections to the second harmonic. 

Advanced ground-based detectors such as Advanced LIGO and EGO will be more 
sensitive at low frequencies. As a result, they will (i) be able to detect stellar mass 
systems at much larger distances; and (ii) be sensitive to signals from more massive 
inspiraling binaries such as very asymmetric, intermediate-mass systems with total 
mass 50 Mq ^ M < 1000 Mq. In the time domain it is natural to terminate waveforms 
at a time determined by F{t) = /lso, where /lso is the orbital frequency at the 
last stable orbit. In the SPA this will roughly correspond to cutting off the kth 
harmonic at / = kfiso- Thus, higher harmonics might still enter the bandwidth 
even if the leading-order, k = 2, harmonic does not, and the resulting SNRs can be 
considerable. The consequences are twofold. At relatively small distances (~ 100 
Mpc), the detector's mass reach will increase significantly, by about a factor of two 
in Advanced LIGO and a factor of three in EGO. This would allow for the detection 
of intermediate-mass inspirals with much higher total mass than with the restricted 
waveforms. Secondly, over cosmological distances, restricted waveforms underestimate 
detection rates of intermediate-mass inspirals by at least a factor of 20. 



Our results can be extended in various directions. 



• Comparison with other waveforms. We have investigated how simulations involving 
amplitude-corrected waveforms both for templates and simulated signals would com- 
pare with simulations involving only restricted templates and "signals". Since the 
(2.5, 3.5)PN waveforms are presumably the closest approximations we have to real 
signals, it would be interesting to see how effective the standard Taylor, EOB, Pade, 
and BCV waveforms are in detecting them. 



15 



Parameter estimation. Some consequences of other harmonics for parameter esti 



mation were studied by HeUings and Moore for the case of LISA. A systematic 
treatment indicating quahtative differences in parameter estimation with restricted 
and non- restricted waveforms for ground-based detectors will appear in [l^, and a 
similar study for LISA is in preparation. 

Testing general relativity. The post-Newtonian phasing of alternative theories of grav- 
ity (in particular scalar-tensor and massive graviton theories) has been worked out 
to IPN order by Berti, Buonanno and Will |38|, who also studied the accuracy with 
which these could be distinguished in the context of LISA. Another recent proposal 
by Arun et al. [s^ exploits the interdependence of the parameters ipi in the phasing of 
waveforms (see Appendix A) to look for deviations from general relativity. It would 
be of great interest to devise similar tests that also take amplitude corrections into 
account. 
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Appendix A: Structure of the amplitude-corrected waveforms in the stationary 
phase approximation 



The expression for the phase is (12.91) . where the coefficients ipi are given by^ 0,111, [13] 



^0 
^5 
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20 
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254016 
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-V 



21 

74045 2 

V 

756 ' 



(5.1) 



with 7 = 0.5772 . . . the Euler-Mascheroni constant. The parameters A and 6 were unknown 
until recently because of regularization ambiguities; they have been determined to be A = 
~loio ^ ~ 0' completing the phasing formula up to 3.5PN. 

In the expressions (12. 8p for the h^''^ we take the "frequency sweep" F to be 



48 



743 11 



336 



7] ] (27rMF)2/3 + 47r(27rMF) 



/34103 13661 



+ 



59 



18144 2016 



V + —V {2ttMF) 



4/3 _ 



/41597r 1897r 
V 672 ^ ~ 



-V 



) {2TiMFfl% 



(5.2) 



We remind the reader that our conventions are shghtly different from what one often finds in other papers, 
as explained in Section |TT1 Frequency dependence in equations (|2.9p . (|5.f p . and (|5.2p is to be read with 
this caveat in mind. 
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where Ai = Mrf/^ is the chirp mass. To 2.5PN order, 



mf)) 

^TT A/i 

48 (27rA</)iV6 [1 + Sii2nMff/' + 53/2(2vrM/) + 52(27rM/)^/3 + S,/2{2itM ff^''] 

(5.3) 



where 

1 /743 11 

S3/2 = -27r, 

7266251 18913 1379 , 

>J2 = 1 V H V , 

8128512 16128 ' 1152 ' ' 

(5.4) 

The coefficients H^+^^ in ([211]) can be found in [J] , and we will not repeat them here. They are 
linear combinations of sines and cosines of multiples of the orbital phase The prefactors 
in these expressions will be denoted C^^'^^ and -D^x\ the former being the prefactors of 
cos(n\E') and the latter the prefactors of sin(n\E^). Now define^ 



P(„,,) = sign (Cf + dr^F^ J (Cf + d^'-^rX + + Di-'^F," ' 



(5.5) 



and 



^in,s) = tan M --J—— . (5.6) 



X 



The harmonics contributing to the waveform up to 2.5PN have the following structure. 

/.«(/) = :^yjvr-^/^(2/)"^/^[e-^v.a^v.)P(,,/,)(2vrM/)V3 

+ (e-*'^(^-^/^)P(i,3/2) + e-^^(i.V2)P(,,i/2)5i)(27rM/) 
+ (e-*'^(i'^)P(i,2) + e-^^(M/2)P(^^^/2)53/2)(27rM/)^/3 

+ (e-*'^(^'^/^)P(i,5/2) + e-^^(i'«/^)P(i,3/2)5i + e-*'^(M/2)P(^^^/2)52)(27rM/)5/3] 

X ^(/l50 - /) exp [z (27r/t, - 7r/4 + V^(/))] , (5.7) 



Up to and including 1.5PN order, the h'^^ involve cosines only and the H^^"^ involve sines only, so that 



D^"'"'' = and C'^-'' = until 2PN order. 
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TT 



-2/3 |'g-j¥'(2,0) 
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r V 48 

+ (e-^^(^>i)P(2,i) + e-*'^(^-")P(2,o)5i)(7rM/)2/3 

+ (e-'^(^>3/^)P(2,3/2) + e-*'^(^'«)P(2,o)53/2)(vrM/) 

+ (e-i^(2,2)p^2,2) + e-^'^(2-i)P(2,i)5i + e-'^(2>o)P(2,o)^2)(vrM/)^/3 

+ (e-*^(^''V^)P(2,5/2) + e-^'^(^-^/^)P(2,3/2)5i + e-*'^(^^^)P(2,i)^3/2 + e-^'^(^'0)P(2,o)^5/2)(vrM/)5/3] 
X 9{2fLso - /) exp [z (27r/t, - 7r/4 + 2^(//2))] , (5.8) 



/,(3)(/) = 3-1/2 




-2/3, 



4g.. (2//3)-^/6 [e-^^(3.V2)P(3^,/2)(2vrM//3)i/3 

+ (e-'^(^>^/^)P(3,3/2) + e-*'^(3,i/2)P(3^^/2)5i)(27rM//3) 

+ (e-M3,2)p^3^2) + e-^^(3.V2)P(3^^/2)^3/2)(2vrM//3)^/3 

+ (e-'^(^'^/^)P(3,3/2)5i + e-'^(3>v^)P(3,i/2)^2)(27rM//3)^/='] 

X 9{3hso - f) exp[z(27r/t, - n/4 + 3^A(//3))], (5.9) 

/^^^H/) = 4-i/2:^yi7r-2/3(//2)-V6[e--(M)P(,^,)(,,M//2)2/3 

+ (e-*^(^'^)P(4,2) + e-*'^(4-^)P(4,i)5i)(7rM//2)^/3 

+ (e-*^(^'^/^)P(4,5/2) + e-*^(^.^)P(4,i)53/2)(vrM//2)5/3] 

X ^(4^50 - /) exp[z(27r/t, - 7r/4 + 4^(//4))], (5.10) 

h^'\f) = 5-i/2^^-7r-2/3(2//5)-^/^[e-^^(^.3/2)p^^^3/^^(2;rM//5) 

+ (e-*'^(^'^/^)P(5,5/2) + e-^'^(-^/2)P(5,3/2)5i)(2vrM//5)^/3] 

X ^(5/l50 - /) exp[z(27r/te - 7r/4 + 5V^(//5))], (5.11) 



/,(6)(/) = 6-1/2 



A^5/6 




7r-2/=^ (//3) -7/6e-'^(M) P(g^2) (vrM//3) 



r V 48 

X e{6hso - f) exp[z(27r/t, - 7r/4 + 6^/;(//6))], 



(5.12) 



/,(7)(/) = 7-i/2:^yj^"2/3(2//7)-7/6g-i^(,,.)p^^^^^^27rM//7)^/=^ 
X 0{7fLso - f) exp[z(27r/t, - n/4 + 7V;(//7))]. 
These can in turn be written as 



r V 48 



X exp 



J]C(fcW2)(2vrM/) 
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(5.14) 
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In the main text we need the C(fc,m/2) for a (1,3.5)PN waveform. The non-zero coefficients 
appearing in such a waveform are given by (with c = cos(t) and s = sin(i)): 



C(i,i/2) = 2-^/'5e-^^(i>i/^)P( 



(1,1/2) 



2-7/6sign --^1 -4r/ (5 + c2)F+ 
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64 16 



X exp 



-i tan" 



-1 



-6cFx 



(5 + c2)F+;_ 



(5.15) 
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C(2,l) = 2-1/22-2/3 ^g-i^(,,i)p^^^^^^g-i^(,.0)p^^^^^ 

= 2-^/6^sign [(19 + - 2c^ - (19 - llc^ - 6c^)r/) F+] 
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Appendix B: Proof of the inequality (|3.6p 



From the expressions 05. 151) . 05.161) . 05.17p . and 05.181) for C(i,i/2), C(2,o), C(2,i), and C(3,i/2), 
one has 

|C(1,1/2)P + IC(3,l/2)P + 2Re[C(\o)C(2,l)] 
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'743 11 
336 ^ ~i 

where we have used that 

sign [-(1 + c2)F+] sign [(l9 + 9c^ - 2c^ - (19 - llc^ - 6c^)r/) F+] = -1 

Applying basic trigonometric identities in the third term of (15.201) . 

|C(i,i/2)r + |C(3,i/2)r + 2Re[C(2,o)C(2,i)] 

9 



(5.20) 



(5.21) 
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The RHS can be written as a linear combination of and , and sufficient conditions for 
fl5.22p to be negative definite are tliat the coefficients in this combination are negative for 
any c G [—1, 1] and t] G [0, 1/4]. It is straightforward but tedious to show that this is the 
case. Note the presence of the prefactor (1 — Arj) in several of the positive contributions to 
these coefficients; this explains the strong //-dependence of the SNR reduction. 

The origin of the inequality fl3.6l) . and similar inequalities one would encounter at higher 
order, can be understood as follows. The cross-term in the LHS of f l5.22p is of the form 
-^^[Cffc (m-i)/2)C(fc,(m+i)/2)], which has a large (in absolute value) contribution that looks like 

yV Slgn[G4. JSlgn[G4. \ \^{k,{m-l)/2)^{k,{m+l)/2)\ 

X COs(v9(fc,(m-l)/2) - V^(fc,(m+l)/2))), (5.23) 

where we have used the notation of Appendix A. Here Af is some positive numerical pref- 
actor. It can be checked that the angles v^{A:,{m±i)/2) always have the same sign and, being 
defined through tan~^, they are in the interval [— 7r/2, vr/2]; hence the cosine will be positive. 
However, C^'^"^ and c^'^^'^^^^'^^ have opposite signs, so that the expression (15.231) is 
negative definite. 



Appendix C: Specifications and projected noise power spectral density for EGO 

The EGO detector is not yet on the drawing boards. Rather, its noise power spectral 
density should be viewed as a summary of what is projected to be possible with steady 
advances in interferometer technology over the next decade or so. For concreteness, a 3 
km arm length was assumed, to be placed underground. The optics are kept at cryogenic 
temperatures (~ 5 K). The loss angle of the mirror substrate and the suspension material 
is approximately 10~^. The injected laser power is 100 W, and the recycling factor of the 
recycling mirror is 50; the finesse of the main cavities is about 600. The radiation pressure 
noise is suppressed using a scheme suggested by Courty et al. 13], in which a control cavity 
is coupled to each mirror, storing a lower light power. The seismic filtering is realized by 
combining a passive system similar to VIRGO's "super-attenuator" with an active system as 
in Advanced LIGO. Above frequencies of a few Hertz the resulting PSD is well-approximated 
by 



Sh{f) = So 



1 + hix + h2x'^ + b^x^ + h^x^ + h^x^ + h^x^ 

1 + CiX + 02^2 -|- C^X^ + CiX'^ 



(5.24) 



where x = f / fo with /o = 200 Hz, and 5*0 = 1.61 x 10"^^ Hz^^ One has 



bi = 31.18, b2 

Cl 



= -4.05, 
ai = 185.62, 
-64.72, 63 = 52.24, 



13.58, C2 



-36.46, 



P2 = -0.69, 
a2 = 232.56, 

64 = -42.16, 65 = 10.17, be 
C3 = 18.56, C4 = 27.43. 



11.53, 
(5.25) 



The strain sensitivities for Advanced LIGO and EGO are plotted in Fig. HI 
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Finally, the lower cut-off frequencies fs for EGO as well as Initial and Advanced LIGO 
are chosen in such a way that for any (untruncated) restricted PN waveform ho, 



'j;^\ho{f)M/s,{f) 

r\ho{fWdf/S,{f) 



1/2 



j;;r"'df/s,{f) 

J^^f-VMf/S,{f) 



1/2 



> 0.99. 



(5.26) 



Picking the largest value of fs satisfying this requirement and rounding off yields the usual 
/, = 40 Hz for Initial LIGO and /, = 20 Hz for Advanced LIGO, while = 10 Hz for EGO. 



■Advanced LIGO 
EGO 




Figure 4: The strain sensitivities of Advanced LIGO and EGO as functions of frequency. 
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